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1. Introduction 

Conformal or near conformal behavior in field theory lie at the heart of many challenging 
theoretical and phenomenological problems. For example models that seek to replace the Higgs 
mechanism of the Standard Model with a new strong gauge dynamics at the TeV scale often invoke 
conformality as an explanation of the flavor hierarchies. While lattice gauge theory in principle 
provides a useful approach to explore such non-perturbative dynamics, conventional lattice meth- 
ods for theories that are parametrically close to conformal theories are difficult precisely because 
of the growing separation of the length scales between the UV and IR. Here we explore a new 
technique that replaces the traditional Euclidean lattice in favor of one suited to Radial Quantiza- 
tion. Radial quantization has a long history starting with the observation that the early covariant 
quantization of the 2-d conformal string action was given as a radial quantized system with the 
Virasoro Lq operator replacing the Hamiltonian. In 1979 Fubini, Hanson and Jackiw [|l|] suggested 
radial quantization of field theory in higher dimensions and later in 1985 Cardy suggested lattice 
implementations in general dimensions 

For an exactly conformal field theory, the idea is straight forward. The flat metric for any 
Euclidean field theory on MP can obviously be expressed in radial co-ordinates, 

d^s = dx^dx^' =ry{dt^ + dQ.l-\) , (1-1) 

where t = log(r/ro), introducing an arbitrary reference scale ro, and where dQ.jj_i is the metric on 
the S^^' sphere of unit radius. However in the case of an exactly conformal field theory, a local 
Weyl transformation will also remove the conformal factor, exp[2f], from the Lagrangian of the 
quantum theory. The resultant theory is mapped from the Euclidean space to a D dimensional 
cylinder, M x S^^^ A simple intuitive illustration of this map begins with the exact two point 
function for a primary operator with dimension A, 

((/)(^l)0(.T2))= I (1.2) 

\Xi — X2\ 

and then converts it to radial form, 

as t = log(r2) — log(ri) oo. The factors on the left are the Weyl factors absorbed into the field 
redefinition of operators for radial quantization. The angular dependence projected onto spherical 
harmonics give rise to integer spaced descendants: A/ = A + /. 

Our goal is to develop numerical methods to solve conformal quantum field theories as the 
infinite refinement limit of a lattice regularization on M x S^^^ If the action is real, one can 
solve the latter numerically by Monte Carlo methods, extract quantitative features and test the 
mathematical question of convergence to a universal continuum limit. If this is possible, a potential 
advantage is that a lattice with T sites mt = log(r/ro) represents an exponential scale separation 
as function of T relative to conventional Euclidean finite lattice on R^. 
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2. Lattice Implementation 

As a test of this idea, we present first results for the 3D Ising model at the Wilson-Fisher critical 
point. The largest discrete subgroup of the isometrics of M x are achieved by a uniform lattice for 
the non-compact M co-ordinate and an icosahedral approximation to the compact sphere, S^. The 
icosahedron has 12 vertices and 20 faces given by identical flat equilateral triangles as illustrated in 
Fig. [I[ Its symmetry group //, is a 120 element subgroup of 0(3). The angular momenta / = 0, 1,2 
representations of 0(3) remain irreducible representations under //,. There is enough symmetry 
to isolate a scalar primary state and at least two of its immediate descendant states. To refine the 
lattice on the icosahedron each face is subdivided into equilateral triangles. The sites on each 
icosahedral surface at fixed t are connected to the corresponding sites on the neighboring surfaces 
at f ± 1. 

The partition function is of the usual form, 

2 = rr^gIf,.v^f^(^^)f^(f+l'-^)+Ir,(.v>.>^C^(^-^)<7(f,3') ^ (2.1) 

where (xy) denotes a nearest-neighbor pairs on each icosahedral shell and f = 0, • • • T — 1 sums over 
the radial co-ordinate. The trace is the sum over the Ising spin, o{t,x) = ±1, on each site {t,x). 
For finite s, the logarithm of the transfer matrix along the cylinder is a regularized representation of 
the dilatation operator. To get information on the spectrum of the transfer matrix, it is convenient 
to compactify the infinite axis of the cylinder to a circle with periodic boundary conditions on the 
spins. 




Figure 1: On the left is the s = \ icosahedral approximation to the sphere and on the right the 5 = 8 
equilateral triangle refinement of the icosahedron , illustrated by projecting each vertex at fixed angles from 
the center of the icosahedron to the unit sphere. 

To approach the Wilson-Fisher conformal field theory in the continuum limit, we need to tune 
j3 to the critical point. The relative scale between the longitudinal and transverse lattice ("speed 
of light" ) is fixed by the integer spacing of descendants of the primary operators. There are no 
other free pai^ameters. For example the leading primary operator odd under Z2 has a sequence of 
descendants A/ = l/2 + T]/2 + /for/ = 0,l,.... The first 3 states were clearly identified with modest 
calculations, verifying the integer spacing and determining the anomalous contribution T]. 
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3. Numerical Results 

The number of sites on one icosahedral shell is 2 + lO^^. We chose our cylinders to have 
lengths which scale with the refinement T = ps. To locate the critical point, we used aspect ratios 
p = 4, 8 while for computation of the magnetization correlation functions we used only p = 8. 

The critical point jS^ was determined first by constructing a sequence of pseudo-critical j3- 
values defined as matching points of the Binder cumulants[||], 



UiP,s,p) 



1 



3(M2)2 



(3.1) 



for s 



rsp at consecutive p-values. This was done for several values of r in the range 1.5 < r < 
5. We also numerically obtained the subleading terms in the approach to the fixed point predicted 
by the Renormalization Group. Subsequently we improved the estimate by performing a global fit 
of the scaling relation, 



U{l5,s,p)^U{P,,oo,p) + aiip)[P-P,]s'^'+bi{p)s 



(3.2) 



to many independent simulations (p = 4,8), a subset of which appear in Fig. ^ of the scaling 
relation where the exponents v and ft) were held constrained to agree with the published values [Q]. 
We find pc = 0.16098698(2), U{pc,^A) = 0.3040(2) and U{l3c,°°,8) = 0.1876(2). 
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Figure 2: Determining j3c from the Binder cumulants J/(s) = 1 — (M^)/3(M^)^ near the pseudo-critical 
point for two different aspect ratios p and increasing values of s. 



We started our study employing the Swendsen-Wang cluster algorithm but switched to the 
more efficient single cluster Wolff algorithm [^. For our final results on the spin-spin con^elation 
function, we generated ensembles at j3 = 0.160987. Each independent run was thermalized using 
2048 sweeps of the Wolff algorithm followed by 8192 sweeps with one estimate of the spin-spin 
coiTclation function after each sweep. We defined a sweep to be 19^/2 Wolff cluster updates which 
leads to the average number of spins flipped each sweep equal to the total volume. All results for 
a given run are then averaged together to form a single blocked estimate and many thousands of 
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independent block estimates are combined to form the ensemble at each s. The jackknife method 
was used to estimate eiTors. 

We projected the spin-spin correlations function on spherical harmonics, 

Q{t)= £ YUa,){a{t + to,x)o%,y))YUQ.y) (3.3) 

where Yim{Q.x) is spherical harmonic evaluated at the angular position of the site x, weighted by 
1/3 of the area of the adjacent spherical triangle projected on the unit sphere as illustrated in Fig. |l]. 
This represents a finite element approach giving an improved approximation to orthonormality of 
the discrete spherical harmonic. As we are only interested in the rotationally-invariant part of the 
correlation function on any given lattice, we have summed over the azimuthal quantum number, m. 

In addition we found it very useful to evaluate the spin-spin correlation function using the 
momentum space single cluster improved estimator method |^]. The connected correlator for the 
lowest mass discrete eigenstate of transfer matrix on our periodic lattice is exactly represented by 
a single hyperbolic cosine, 

Q(0=A,cosh(-At/(?-r/2)), (3.4) 
at discrete values ? = 0, ■ • • , T — 1. We transform this to momentum space, 

r 1 v' itk riA ^ ^ (1 - e'^'^ sinh(At/) 

J smh^(/X//2) + sm^(^/2) 

where k = 2nq/T with <7 = 0, • ■ • , T — 1 is the momentum conjugate tot = log r along the cylinder. 
Since our value 01^-^ j3c is slightly larger that the pseudo-critical coupling at any finite s, we expect 
that our / = correlation function will have a small disconnected contribution. This contributes a 
non-analytic term, coSkfl. In momentum space the disconnected piece can in principle be isolated 
by subtracting a smooth extrapolation of Co(^) from k^Otok = 0. 

We found that our data require parameterizing the ground state plus at least three higher mass 
states to get excellent fits with /Aof < 1 and estimates of ground state masses which are essen- 
tially free of higher state contamination. These represent possible high states either propagating 
forward in the Z2 even or propagation backward in Z2 odd sectors. While the dimension of the 
higher state coming from the Zj even sector might be lower, its mixing will also be proportional to 
the 3 point coupling of the energy operator and two spin operators. Indeed with higher statistics, 
we believe quantitative determination of the higher spectrum is well within reach of this method. 

Once we have determined the jLt/'s, we relate them to the eigenvalues of the dilatation operator 
up to a single unknown constant: /X/ = A^' [Ag + /] where A^' ~ ci/5 as 5 — )• 00. Numerically, we 
find c\ 1.51(1) with the uncertainty dominated by systematic enm. Clearly, we see evidence for 
sub-leading contributions &{\/s^) as well in the left figure of Fig. ^. We then test for the equal 
spacing rule of descendants by examining the ratios, (/I/+2 shown on the 

right in Fig. ^ Using this confirmation, we are able to estimate numerically the scaling dimension 
of the primary operator using ratios. 



pi-^i' l-V 



(3.6) 



as shown in Fig. ^ 
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FigureS: The left figure shows the scaling of A relative to i. The extrapolated value is roughly 1.51(l)with 
the uncertainty dominated by the systematic difference between the two estimates. The right figure tests 
the hypothesis of integer level spacing between primaries and descendents by comparing the level spacing 
between the first and second descendents to the spacing between the primary and its first descendent. We fit 
to a linear function and find an intercept of 0.994(4) and slope of 0.0(2) with j^/dof = 0.43 for 1 1 dof. 



We end up with what we feel is a quite conservative estimate of t] = 0.034(10), consistent with 
other estimates [Q]. Moreover this method can be extended to include additional primary operators 
in both the Z2-0AA and Z2-even sectors as well as direct test of the restoration of full conformal 
symmetry for 2- and 3-point correlators. 



4. Discussion 

We have presented a simple example of lattice radial quantization for the 3D Ising model. This 
raises important questions and suggests further applications. 

For the 3D Ising radial quantization, we have chosen a simple approximation to the sphere by 
a triangular refinement of the icosahedron with "flat" sides. This is reflected in the action by assign- 
ing equal weights on all the nearest neighbor links. Viewed in the language of Regge calculus, this 
geometry has all the curvature concentrated at 12 exceptional vertices of the underlying icosahe- 
dron that are bounded by 5 rather than 6 triangles. We are implicitly making the conjecture that by 
virtue of maintaining exact icosahedral symmetry the continuum limit gives back the symmetry of 
S^, indeed the full conformal group. Our modest numerical results to date support this conjecture, 
but we are undertaking more stringent numerical test on the spectrum and correlators. It is an open 
question whether the conical singularities at the vertices of the icosahedron are iiTclevant to the 
continuum spectrum. If necessary one might introduce improved triangulation of the metric on S^, 
similar to our improvement of the y/„, weights in our correlation measurement. On the other hand, 
it is also interesting to ask if the simpler geometry of a refined cube is adequate, since applications 
to 4D gauge theories are probably easier to formulate on concentric 3D hypercubes. 
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Figure 4: The scaling exponent of lowest Z2-0AA primary operator vs. \/s. The lower (red) points using 
the primary and first descendant have disconnected contributions that have not been fully determined. A 
linear extrapolation from upper (blue) points, using the first two descendants which have no disconnected 
contributions, gives Aq = 0.517(5), is consistent with the best published estimate 0.5182(3) 

The next simplest model beyond the 3D Ising model to consider is the 3D 0(N) model, which 
because of the analytical results in the large N limit is an excellent test-bed for the method. We 
are considering generalization to include gauge fields and fermions in 3D and 4D. Each of these 
steps will require careful consideration to make sure that there are no obstructions to taking the 
continuum limit. There maybe subtle issues for example with fermions on a spherical manifolds 
and potentially relevant conformal symmetry breaking operators the don't vanish in the continuum. 
Even more challenging are theories that are not quite conformal where the dilatation operator is no 
longer a conserved quantity, such as those exhibiting asymptotic freedom in the UV and softly 
broken conformality near an IR fixed point. 
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